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ABSTRACT 

We report the discovery of an eclipsing companion to NLTT 41135, a nearby M5 dwarf that was 
already known to have a wider, slightly more massive common proper motion companion, NLTT 41136, 
at 2'.'4 separation. Analysis of combined-light and radial velocity curves of the system indicates that 
NLTT 41135B is a 31 — 34 ± 3Mj up brown dwarf (where the range depends on the unknown metallicity 
of the host star) on a circular orbit. The visual M-dwarf pair appears to be physically bound, so the 
system forms a hierarchical triple, with masses approximately in the ratio 8:6:1. The eclipses are 
grazing, preventing an unambiguous measurement of the secondary radius, but follow-up observations 
of the secondary eclipse (e.g. with the James Webb Space Telescope) could permit measurements of 
the surface brightness ratio between the two objects, and thus place constraints on models of brown 
dwarfs. 

Subject headings: binaries: eclipsing - stars: low- mass, brown dwarfs - stars: individual (NLTT 41135) 



1. INTRODUCTION 

The formation mechanism of brown dwarfs, star-like 
objects whose masses are too low to sustain hydrogen fu- 
sion in their interiors, has remained a theoretical puzzle 
ever since they were first detected. They pose a problem 
for standard star formation theories relying on gravita- 
tional collapse of molecular cloud material because their 
masses are substantially smaller than the typical Jeans 
mass in such a cloud. 

Although considerable theoretical and observational 
effort has been expended, there is still no clear con- 
sensus on how brown dwarfs form. Popular theoreti- 
cal ideas can be divided into three major sub-groups. 
The first relies on finding some mechanism to produce 
a sufficiently low-mass core ab initio, while circum- 
venting the Jeans mass problem. This can be done, 
for example, by turbulent fragmentati on of cloud ma- 
terial (e.g. iPadoan fc Nordlundl 120041) . fragmentation 
and s ubseq uent collapse of molecular cloud cores (e.g. 
iBossl 120021) or erosion of core s by ionizing radiation 
from nearby massive stars (e.g. IWhitworth fc Zinneckerl 
2004). The second posits that brown dwarfs form as 
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part of multiple systems (which collapse from a single 
large overdensity, thus circumventing the Jeans mass 
problem), and are then ejected before they can ac- 
crete sufficient m ass to become normal hyd rogen burn- 
ing stars (e.g. IReipurth fc Clarke! 120011). Finally, 
disk fragmentation mechanism s (e.g. iRice et ail 120031 : 
IStamatellos fc Whitworthll2009[ ) are proposed, by which 
stars and brown dwarfs form by gravitational instability. 
Similar mec hanisms hav e bee n proposed for plane t for- 
mation (e.g. IBossll2006D . and IKratter et all (|2010l) note 
that disk fragmentation in the outer regions of proto- 
planetary disks naturally leads to the formation of brown 
dwarfs. 

Multiple systems offer an important clue as to the na- 
ture of the star formation process and could give an in- 
sight into the problem of brown dwarf formation. Ob- 
servations indicate that while multiple systems are ex- 
tremely common for solar type stars, with approximately 
65% of systems found to harbor o ne or more components 
(e.g. iDuquennov fc Mavorlll99i|) . the multiplicity frac- 
tion appears to decline through M spectral types to the 
brown dwarf domain, with only 42% of field M-dwa rf sys- 
tems being multiple (e.g. iFischer fc Marcv| [T992). The 
fraction could be as low as 10% for late-M and L-dwarfs, 
although there is some disagreement in the literature re- 
garding corrections for the effects of su rvey biases, par- 
ticularly in binary semimajor axis (see [Burgassc r et al.l 
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120071 and references therein for a review). One forma- 
tion mechanism that naturally generates a low multiplic- 
ity fraction for br own dwarfs is the ejectio n hypothesis 
mentioned above (Rcipurth fc Clarke! 120011 ) since it re- 
lies on dynamical interactions, which would likely dis- 
rupt most binary brown dwarf systems. However, this 
mechanism has difficulty explaining other observational 
evidence, s uch as the existence of di sks around brown 
dwarfs (e.g. Uavawardhana et al.ll2003l ). which would also 
be disrupted by dynamical interactions. 

In addition, binary systems have long been a corner- 
stone in our understanding of the fundamental properties 
of stars because gravity provides one of the few model- 
independent means to precisely measure the mass of a 
star for comparison with other fundamental properties, 
thus allowing us to test stellar evolution theory. Double- 
lined eclipsing binaries (EBs) are particularly useful in 
this regard, where combined analysis of light curves and 
radial velocity information can determine masses and 
radii for both stars to very high p recision (iTorres et al.l 
2010; see also the classic review by IAndersenlll99lD . 

Results of recent dynamical mass measurements for 
brown dwarf binaries indicate that the there may be sub- 
stantial discrepancies between the predictions of evolu- 
tionary models and the observed masses and luminosities 
of the objects, by factors of up to 2 — 3. There has been 
some debate in the literature as to the nature of these 
discrepancies, with the earliest example findin g the lumi- 
nosit ies were overpredicted for a given mass (jClose et al.1 
120051 : but this may have been du e to an error in the 
spectral type, see lClose et al.ll2007D , whereas other stud - 
ies found they are underpredicted (IDupuv et "all [ 2009). 
or indicate no dis crepancy (e.g. iLiu et al.l 120081 : and 
IDupuv et all 120091 for an intriguing object at the L/T 
transition) . By making measur ements of an unpre c edent - 
edly large sample of objects, iKonopackv et al.1 ()2010l ) 
conclude that the discrepancies are a function of spectral 
type, with late-M to mid-L systems having luminosities 
overpredicted by the theoretical models, and one T-type 
system having a luminosity which is underpredicted, in- 
dicating that the sense of the discrepancies may reverse 
moving to later spectral types. Furthermore, the evo- 
lutionary and atmospheric models may indeed be in- 
consistent with one another (e.g. IDupuv et all 12009: 
IKonopackv et al.ll20 10t), indicating systematic errors in 
one or both classes of models. 

The only well-characterized double lined brown dwarf 
EB is perhaps even more puzzling, showing a reversal 
in effective temperatures of the two obj ects, with the 
more mass ive brown dwarf being cooler (|Stassun et al.l 
120061 12007T) . This is thought to be due to magnetic ac- 
tivity in this extrem ely young (~ 1 — 2 Myr) system 
(jChabrier et al.l 120071 ). What is clear is that even more 
dynamical measurements of brown dwarfs are needed, to 
place these results on a firmer footing and to probe a 
wider range of parameter space. 

In this work we present the detection of an eclipsing 
brown dwarf orbiting the lower-mass member of a known 
field M-dwarf visual double, NLTT 41136 and 41135, 
which we show are highly likely to be physically asso- 
ciated. The observation of eclipses constrains the orbital 
inclination, and thus allows the secondary mass to be as- 
certained through radial velocity measurements with es- 
sentially no sin i ambiguity. The properties of the parent 



M-dwarf can also be constrained from the eclipse mea- 
surements. Such an object represents an opportunity to 
constrain formation and evolutionary models of both M- 
dwarfs and brown dwarfs. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. MEarth photometry 

Eclipses in NLTT 41135 were initially detected from 
data taken during 2009 February to May (inclusive) 
as part of routine operation of the MEarth observa- 
tory, a system designed primarily to search for transit- 
ing super-Earth exoplanets orbiting around the nearest 
2000 mid to late M-dwarfs in the northern hemis phere 
(|Nutzman fc Charbonneaull2008l; llrwin et al.ll2009D . Ex- 
posure times on each field observed by MEarth are tai- 
lored to achieve sensitivity to the same planet size for the 
assumed parameters of each target star, and were 151s 
for the field containing NLTT 41135. 

Data were reduced using the standard MEarth reduc- 
tion pipeline, which is at present ne arly identica l to th e 
Monitor project pipeline described in llrwin "et^lpOOl . 
We used an aperture radius of 10 pixels (7'.'6); note that 
the MEarth telescopes were operated slightly out of fo- 
cus in the 2008 September to 2009 July observing season 
to increase the number of photons that can be gathered 
before saturation of the detector for our brightest targets. 

This star was identified as a high-priority candidate 
from a box-fitting least squares (BLS) transit searc h (fol- 
lowing closely the methods of iBurke et al.l [2006;) , per- 
formed on de-trended, me dian filtered light curve s. We 
used a 3-day median filter (Ai grain fc irwinl l2b04') to re- 
move stellar va riability, followed by a trend filtering al- 
gorithm (TFA; iKovacs et al.l 120051 ) based on searching 
for optimal linear combinations of comparison star light 
curves to best remove any systematic effects in the tar- 
get light curve. The photometry indicated a period of 
approximately 2.8 days, and a 2% eclipse depth, with no 
secondary eclipses visible in the MEarth data. We im- 
mediately switched to a follow-up mode after the initial 
detection, observing the object at the highest possible 
cadence (approximately 3 minutes including overheads) 
during the night of UT 2009 May 25 when an eclipse was 
predicted to occur. This was confirmed with high signif- 
icance, although any measurements of the eclipse shape 
were severely corrupted by the target crossing the merid- 
ian during the eclipse, which necessitates a temporary 
halt in the observations and rotation of the telescope op- 
tics and detector system through 180° relative to the sky 
with our German Equatorial Mounts. In practice, the 
telescope guiding takes some time to recover after this, 
and photometry is corrupted during the recovery period. 
Therefore while we were able to confirm the reality of the 
event, the MEarth data are not useful in determining the 
eclipse shape or physical properties of the object. 

Due to the close proximity of NLTT 41135 and 41136, 
these stars are unresolved in the MEarth observations, 
and indeed are also unresolved or improperly resolved 
in many well-known literature sources for measurements 
of bright stars. It was initially not clear which mem- 
ber of the M-dwarf pair was undergoing eclipses, so we 
pursued photometric (resolved light curves) and spectro- 
scopic (radial velocity) methods to determine the identity 
of the eclipsing star. The spectroscopic method gave the 
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earliest indication that NLTT 41135 was responsible for 
the eclipses. 

The high proper motion of the pair allows us to con- 
strain the contribution of any additional background 
stars in the photometric aperture that are not co-moving 
with NLTT 41135 and 41136, using previous epochs of 
imaging. We show in Figure [1] a series of three images 
centered on the position of the photometric aperture used 
in the MEarth images. These show that there is a fainter 
star located approximately 7 arcsec to the east and 1 arc- 
sec to the north of the blended pair at the MEarth epoch. 
Using the APM on-line sky catalog, this star is 4 mag- 
nitudes fainter than the blended image of NLTT 41135 
and 41136 in the POSS-1 red (E) plate, and 2.3 magni- 
tudes fainter in the blue (O) plate. While this star is at 
the edge of the MEarth photometric aperture and thus 
will affect the eclipse depths measured from these data, 
we have ensured that this star was excluded from the 
apertures used to analyze the follow-up data described 
in the subsequent sections, and its fainter magnitude and 
bluer color should ensure any contamination in the far- 
red bandpasses used in this work is negligible. 

Our MEarth light curves show variations out of eclipse. 
These have a peak-to-peak amplitude of ~ 0.04 mag, and 
rms scatter of 0.012 mag. There is no clear period, and 
it is not clear which member of the pair is responsible 
for the variations. Our photometry is corrupted some- 
what by the partially-resolved nature of the pair, so we 
regard this measurement as an upper limit to the actual 
variability level, noting that the absence of a clear pe- 
riod may be indicative of some or all of the measured 
variation being spurious. 

2.2. UH 2.2m z-band resolved follow-up photometry 

In order to determine the flux ratio of the M-dwarf 
pair, and hence the intrinsic flux decrement during 
eclipse, we obtained a single, resolved image of the 
system using the Ort hogonal Para l lel T r ansfer Imag- 
ing C amera (OPTIC; iHowell et all l200l iTonrv et al.l 
2005) on the University of Hawaii 2.2m telescope, dur- 
ing the night of UT 2009 June 9. This instrument 
uses two 2048 x 4 096 "orthogonal transfer" CCDs (e.g. 
ITonrv et al.l 119971 ) which allow parallel transfers of the 
charge through the image area in both the x and y direc- 
tions. The pixel scale is O'.'14/pix on the 2.2m telescope, 
giving a field of view of 9.'3 x 9f3. An exposure time 
of 10 s was used in good seeing conditions through the 
z filter. The Heliocentric Julian Date of mid-exposure 
was 2454991.910257, corresponding to an orbital phase 
of 0.23 (close to quadrature) for the NLTT 41135 system 
(see gZHJ). 

This single image was bias-subtracted using the over- 
scan region, and flat-fielded before using our standard 
source detection software to obtain positions and aper- 
ture photometry for both stars in the pair. We mea- 
sured an angular separation of 2'.'40, and a flux ratio of 
^4ii36,z/^4ii35,2 = L3 — 2.1 ±0.2, where we have assumed 
a conservative error based on multiple measurements of 
the image made using different sized apertures. This 
quantity will later be used in the analysis of unresolved 
light curves of the system to extract the parameters of 
the eclipsing binary, NLTT 41135. In this context it is 

3 http : //www . ast . cam. ac .uk/~apmcat/ 




Fig. 1.— Images of NLTT 41135 and 41136 centered on the po- 
sition as measured from the MEarth data. The circle shows the 
approximate position and size of the 7'.'6 (radius) photometric aper- 
ture used to derive our light curves, and the arrow points to the 
position of the fainter star discussed in the text. Data are from the 
first and second epoch Palomar sky surveys as provided by the Dig- 
itized Sky Survey (DSS; top and center panels), and the MEarth 
master image (bottom panel). The two stars are unresolved, al- 
though elongation of the combined image is visible in the POSS-2 
(center) panel. The approximate epochs of the images are 1954.4 
(POSS-1), 1993.2 (POSS-2), and 2009.1 (MEarth). All three pan- 
els have the same center, scale and alignment on-sky, with north up 
and east to the left, covering approximately 2'. 5 in the horizontal 
direction. 

usually referred to as "third light", and shall hereafter 
be labeled with the symbol L3 for consistency with the 
usual nomenclature in the eclipsing binary literature. 

It is important to consider the possible influence of 
stellar variability on the value of L3. Using the results 
stated in i j2.1[ in a similar passband, the most pessimistic 
assumption is to assign all of the measured variation to 
NLTT 41135, giving a contribution to the fractional un- 
certainty in L3 of 4% (rms). Our assumed uncertainty of 
±0.2 (a fractional uncertainty of 10%) in L3 is already 
sufficient to account for this, and the additional contri- 
bution from variability is probably negligible. 

An attempt to obtain a resolved light curve was made 
on the night of UT 2009 June 17. Exposure times of 
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25 s were used through the z filter. We used the orthog- 
onal parallel transfer capability of the OPTIC detector 
to shift charge during the exposure in an attempt to per- 
form a tip-tilt correction of the incoming wavefronts to 
improve the FWHM of the images. However, due to 
technical problems this instead caused severe corruption 
of the point spread functions, which limits the utility of 
the resulting light curve. In addition, during the eclipse 
the guiding lock was lost which caused a short gap in 
the phase coverage. Nonetheless, we were able to extract 
resolved photometry of the two stars by performing stan- 
dard aperture photometry on the resulting images. An 
aperture radius of 4 pixels (0'.'56) was found to give the 
best compromise between flux losses from the aperture 
and cross-contamination of the point spread functions of 
the two stars. We used NLTT 41136 as the comparison 
star to extract differential photometry of NLTT 41135, 
which has the additional advantage of canceling much 
of the cross-contaminated flux between the two objects 
and hence providing a more accurate measurement of the 
true eclipse depth. This light curve is given in Table Q] 

2.3. FLWO 1.2 m z -band follow-up photometry 

Observations centered around the primary eclipse of 
UT 2009 June 20 were obtained using the KeplerCam 
instrument on the Fred Lawrence Whipple Observatory 
(FLWO) 1.2 m telescope. We used the standard binning 
2x2 readout mode, since the pixel scale of 0'.'34 per 
unbinned pixel significantly oversamples the typical see- 
ing at FLWO. The resulting scale was 0'.'67 per summed 
pixel. We used the z filter and an exposure time of 120 s. 
A total of 86 observations were taken, starting approxi- 
mately 1 hour before first contact and finishing 0.5 hours 
after last contact to sample the out-of-eclipse portions 
of the light curve, thus allowing the eclipse depth to be 
properly measured. 

These photometric data were reduced using the same 
pipeline as described in £j 2 . H The FWHM of the stellar 
images was approximately 2'.'5, so we used an aperture 
radius of 8 binned pixels, corresponding to 5'.'4 on-sky, 
to extract aperture photometry of the combined light of 
NLTT 41135 and NLTT 41136, reproduced in Tabled 
Our attempts at PSF fitting photometry on the images to 
extract individual light curves did not yield useful results. 

2.4. Sloan Digital Sky Survey photometry 

Our target lies within the survey area for the Sloan 
Digital Sky Survey ( SPSS). In Data Release 7 (DR7; 
lAbazaiian et al.ll2009| ). the object is resolved in the SDSS 
images, but appears to be blended with a false galaxy 
detection in the band-merged catalogs, so some of the 
magnitudes may be unreliable. We therefore re-derive 
resolved magnitudes for the two stars from the reduced 
images. 

The relevant fpC "corrected frames" were retrieved 
from the SDSS archive. We used the r-band image for 
source detection, since this had the smallest FWHM of 
~ 2.8 pixels. The positions derived therefrom were then 
used to perform aperture photometry on the g, r and z 
band images, using an aperture radius of 2.5 pixels (l'.'O) 
to avoid overlap between the apertures placed on the two 
components of the visual binary. The i image appears to 
be corrupted for NLTT 41136, so we did not attempt to 



derive photometry from it, and the u band image showed 
very low counts for the M-dwarf pair, as expected, so this 
was also omitted, and is likely of limited utility in any 
case due to a well-documented red leak issue with the 
SDSS u filter affecting measurements of very red targets. 

The final calibrated magnitudes are given in Table 
[H We estimate uncertainties in these measurements by 
comparing our measurements of other stars in the image 
with the SDSS DR7 photometry, and allowing for cross- 
contaminated flux between the apertures placed on our 
targets. The photometric errors in g are large in part due 
to the greater FWHM of the stellar images in this band 
compared to r and z, necessitating a larger aperture cor- 
rection. The orbital phase of these measurements for the 
NLTT 41135 system was approximately 0.34 (see §2T8]), 
and the measured z-band flux ratio is 2.03, which is con- 
sistent with the value reported from the OPTIC data in 

a 

Our measurements in g and z for our target stars are 
systematically brighter than those reported in the SDSS 
DR7 catalog, which we suspect to be due to de-blending 
of the false "galaxy" detection from the stellar images. 
The images of our target do not appear to be saturated 
in the SDSS data. 

2.5. FLWO 1.5 m TRES spectroscopy 

Spectroscopic observations of both stars were obtained 
using the TRES fiber-fed echelle spectrograph on the 
FLWO 1.5 m Tillinghast reflector. We used the medium 
fiber (2'.'3 projected diameter), yielding a resolving power 
ofiJ- 44 000. 

We obtained two epochs on each star in the visual bi- 
nary close to the predicted times of quadratures from 
the photometric ephemeris, in order to search for large- 
amplitude radial velocity variability. One hour exposures 
were obtained in June 2009 over a four day period. 

The radial velocities were extracted using a custom- 
built pipeline to rectify and cross correlate the spectra. 
This pipeline was identical to the one used for the Nordic 
Optical Telescope data and will be described in the next 
section. 

The two radial velocities of NLTT 41135 from TRES 
showed variations in phase with the prediction from the 
photometric ephemeris, giving an initial estimate of the 
RV semi-amplitude of K = 12.9 ± 0.5 km s _1 . We 
cross correlated the spectra of the two stars and sub- 
tracted the orbital motion of the low-mass companion of 
NLTT 41135. If the pair of stars are physically bound, 
this relative velocity should be constant over short (few 
day) timescales and close to zero. We found a veloc- 
ity difference of 860 ± 130 m s _1 over 4 days, which is 
consistent with the pair being physically bound, within 
the errors. These results encouraged us to obtain high 
signal-to-noise high-resolution spectra to characterize the 
RV variations of NLTT 41135. 

We note that NLTT 41135 shows the Ha line in emis- 
sion, whereas NLTT 41136 does not, and instead there 
is a hint of an absorption feature at this wavelength in 
our low signal-to-noise spectra. This is reasonably con- 
sistent with expectation for a field-age M-dwarf system 
given the rapid increase in the observed Ha a ctivity frac- 
tion around M4 to M5 spectral types (e.g. iWest et ahl 
l200l . but it may also be indicative of enhanced activ- 
ity on NLTT 41135 due to tidal effects from the close 
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TABLE 1 

UH 2.2m RESOLVED 2-BAND light curve. 



HJD a 


Differential z 


Error b 




FWHM d 


Airmass 


x° 


y c 










(pix) 




(pix) 


(pix) 


2454999.853406 


-0.0502 


0.0039 


-0.567 


6.36 


1.03609 


482.569 


3423.312 


2454999.854022 


-0.0358 


0.0064 


-0.474 


2.95 


1.03604 


474.552 


3421.937 


2454999.854670 


0.0f92 


0.0043 


-0.371 


6.02 


1.03601 


477.566 


3419.809 


2454999.855283 


0.0048 


0.0037 


-0.088 


5.88 


1.03599 


475.421 


3421.553 


2454999.855896 


0.026f 


0.0040 


-0.103 


6.32 


1.03600 


469.411 


3419.771 



Note. — Table [T] is published in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 

a Heliocentric Julian Date of mid-exposure. All HJD values reported in this paper 
are in the UTC time-system. 

b Estimated using a standard CCD noise model, including contributions from 
Poisson noise in the stellar counts, sky noise, readout noise and errors in the sky 
background estimation. 

c Correction to the frame magnitude zero-point applied by the differential pho- 
tometry procedure. More negative numbers indicate greater losses. Please note 
that this has already been applied to the "differential z n column and is provided 
only for reference (e.g. distinguishing frames with large losses due to cloud). 
d Median FWHM of the stellar images on the frame. The plate scale was O'.'14/pix. 
° x and y pixel coordinates on the CCD image, derived using a standard intensity- 
weighted moments analysis. 



TABLE 2 
FLWO 1.2m z-band light curve. 



HJD a 


Differential z 


Error b 


Am c 


FWHM d 


Airmass 


x a 


y c 










(pix) 




(pix) 


(pix) 


2455002.740052 


-0.0000 


0.0012 


0.013 


3.91 


1.12786 


497.620 


524.552 


2455002.741672 


0.0025 


0.0012 


-0.006 


4.00 


1.12904 


497.732 


524.679 


2455002.743188 


0.0000 


0.0013 


-0.019 


4.20 


1.13026 


497.800 


524.485 


2455002.744751 


0.0013 


0.0012 


0.007 


3.87 


1.13162 


497.716 


524.532 


2455002.746302 


-0.0001 


0.0013 


0.016 


3.85 


1.13307 


497.629 


524.695 



Note. — Table [2] is published in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 

a Heliocentric Julian Date of mid-exposure. All HJD values reported in this 
paper are in the UTC time-system. 

b Estimated using a standard CCD noise model, including contributions from 
Poisson noise in the stellar counts, sky noise, readout noise and errors in the sky 
background estimation. 

c Correction to the frame magnitude zero-point applied by the differential pho- 
tometry procedure. More negative numbers indicate greater losses. Please note 
that this has already been applied to the "differential z" column and is provided 
only for reference (e.g. distinguishing frames with large losses due to cloud). 
d Median FWHM of the stellar images on the frame. The plate scale was 
O'.'67/pix. 

e x and y pixel coordinates on the CCD image, derived using a standard intensity- 
weighted moments analysis. 



companion causing it to rotate more rapidly. 

2.6. NOT/FIES spectroscopy 

Precise radial velocities of NLTT 41135 were obtained 
using the Fibre-fed Echellc Spectrograph (FIES) on the 
2.5m Nordic Optical Telescope (NOT) at La Palma, 
Spain. We used the medium resolution hber (l'/3 pro- 
jected diameter) with a resolving power of R ~ 46 000 
giving a wavelength coverage of ~ 3600 — 7400 A. We 
obtained seven high-resolution, high signal-to-noise spec- 



tra of NLTT 41135 during a seven night run in August 
2009. 

The spectra were rectified and cross correlated using a 
custom built pipeline designed to provide precise radial 
velocities for echelle spectrographs. The procedures are 
described in more detail in Buchhave el al. (2010, in 
preparation) and will be described briefly below. 

In order to effectively remove cosmic rays, each ob- 
servation was split into three separate exposures, en- 
abling us to combine the raw images using median fil- 
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tering, removing virtually all cosmic rays. We use a flat 
field to trace the echelle orders and to correct the pixel 
to pixel variations in CCD response, then extract one- 
dimensiona l spectra usin g the "optimal e xtraction" a l- 
gorithm of iHewett etaLl (| 19851 ) (see also iHornd Il986l ). 
The scattered light in the two-dimensional raw image is 
determined and removed by masking out the signal in 
the echelle orders and fitting the inter-order scattered 
light flux with a two-dimensional polynomial. The FIES 
sky fiber was broken at the time these observations were 
taken, so sky subtraction could not be performed. Un- 
fortunately, this means we are also unable to measure 
equivalent widths, e.g. to characterize the activity of the 
star. 

Thorium argon (ThAr) calibration images were ob- 
tained through the science fiber before and after each 
stellar observation. The two calibration images are com- 
bined to form the basis for the fiducial wavelength cali- 
bration. FIES and TRES are not vacuum spectrographs 
and are only temperature controlled to 1/10 of a degree. 
Consequently, the radial velocity errors are dominated by 
shifts due to pressure, humidity and temperature varia- 
tions. In order to successfully remove these large varia- 
tions (> 1.5 km s" 1 ), it is critical that the ThAr light 
travels along the same optical path as the stellar light 
and thus acts as an effective proxy to remove these vari- 
ations. We have therefore chosen to sandwich the stellar 
exposure with two ThAr frames instead of using the si- 
multaneous ThAr fiber, which may not exactly describe 
the induced shifts in the science fiber and can also lead 
to bleeding of ThAr light into the science spectrum from 
the strong argon lines, especially at redder wavelengths. 
The pairs of ThAr exposures are co-added to improve the 
signal to noise ratio, and centers of the ThAr lines are 
found by fitting a Gaussian function to the line profiles 
and a two-dimensional fifth order Legendre polynomial 
is used to describe the wavelength solution. 

Once the spectra have been extracted, a cross correla- 
tion is performed order by order. The spectra are inter- 
polated to a common oversampled log wavelength scale, 
high and low pass filtered and apodized. The orders are 
cross correlated using a Fast Fourier Transform (FFT) 
and the cross correlation functions (CCFs) for all the 
orders co-added and fit with a Gaussian function to de- 
termine the RV. For the FIES spectroscopy typically 10 
orders (of a total of 78) yielded usable cross-correlations, 
where the majority of the orders we did not use are in 
the blue and had very poor signal to noise ratios. For 
TRES, 8-10 orders were used. Uncertainties were esti- 
mated by fitting Gaussians to the individual orders' cross 
correlation functions, and taking the error in the mean of 
all the orders. The template used in this procedure was 
constructed from the target spectra, initially determining 
the RVs by cross-correlating against the best signal-to- 
noise individual spectrum of the target, and then shifting 
and adding all of the spectra to make a single high signal- 
to-noise stacked template, and repeating the correlations 
against this to produce the final RVs. 

The radial velocity measurements of NLTT 41135 are 
reported in Table G3 

Since our radial velocities are measured with respect to 
the target itself, the zero-point is arbitrary. We make a 
separate estimate of the systemic velocity itself (usually 
called 7 in eclipsing binary studies, and we follow the 



TABLE 3 

Barycentric radial velocity measurements 
of NLTT 41135. 



HJD a 


v b 


<7„ 


Sccing c 




(ms- 1 ) 


(ms^ 1 ) 


(arcsec) 


2455048.392461 


22098.6 


31.5 


1.0 


2455049.390174 


-83.2 


27.4 


0.9 


2455050.385098 


6408.9 


48.8 


0.9 


2455051.384396 


21245.8 


25.4 


0.8 


2455052.383913 


-2062.3 


33.4 


0.7 


2455053.381672 


9464.6 


38.2 


1.1 


2455054.373863 


19808.2 


52.7 


1.8 



a Heliocentric Julian Date of mid- 
exposure. All HJD values reported in this 
paper are in the UTC time-system. 
b Relative radial velocity. The zero-point 
is arbitrary and corresponds to the origi- 
nal reference spectrum (HJD 2455049.39). 
An estimate of the actual 7 velocity is re- 
ported in the text and given in Table [3] 
c Seeing estimates for a wavelength 
of 500 nm and an airmass of 1.0 from 
the Differential Image Motion Mon- 
itor (DIMM) operated by the Isaac 

Newton Group on the same site. See 

http : //www. ing. iac . es/Astronomy/development/seeing/ 
tor more information on this instrument. 

same notation in this work) using the Ha emission lines 
in NLTT 41135. These were fit with Gaussian functions, 
and the resulting spectroscopic orbit was analyzed in the 
same way as described in to determine the 7 velocity 
and its error, reported in Table O 

We do not see evidence for significant rotational broad- 
ening in the observed cross-correlation functions, where 
the formal FWHM of the cross-correlation peak was 
w 12 km s -1 . This implies a FWHM of 8.5 km s" 1 in the 
individual spectra as we used the target as its own cross- 
correlation template. Formally, given the stated resolv- 
ing power, this would imply a contribution of « 5 km s _1 
from the star, although the uncertainty in this value 
is large and we do not believe the small difference in 
FWHM is significant in practice, as other sources can 
inflate the measured cross-correlation function width. 
Unfortunately, we lack observations of suitable slowly- 
rotating M-dwarf templates to make a robust measure- 
ment of the rotational broadening. However, the lack 
of a coherent, periodic modulation in the MEarth pho- 
tometry (see £|2.1[) also indicates probable low rotation, 
where usini = 5 km s" 1 would imply a rotation period 
of < 2.1 days. 

2.7. Summary of the photometric and astrometric 
system properties 

Table 3] summarizes the known system properties, from 
the literature (princip ally the proper motion survey of 
iLepine fc~ ghara 200^ used to select the target stars for 
the MEarth survey) and the SDSS photometry. A source 
appears in the 2MASS point source catalog, but we do 
not report these measurements here since the photome- 
try is flagged as having a very poor point spread function 
fit and may therefore be unreliable. Comparing the pre- 
dicted and measured J-band magnitudes indicates that 
the source detected in 2MASS is most likely NLTT 41136 
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TABLE 4 

Summary of the photometric and 
astrometric properties of nltt 41135 and 
NLTT 41136. 



Parameter 


NLTT 41136 


NLTT 41135 


a2000 a ' b 


15 h 46 m 04 s .41 


15 h 46 m 04 s .26 


5 2 ooo a ' b 


+04°41'31'.'2 


+04°41'30?0 


jX a COS 8° 


0"156 


yr" 1 




-0'.'284 


yr" 1 


9 


16.261 ± 0.06 


17.422 ± 0.06 


r 


14.927 ±0.03 


16.036 ± 0.03 


z 


12.374 ± 0.03 


13.145 ±0.03 


Spectral type 


M4.2±0.5 


M5.1 ± 0.5 


T cff c 


3340 ± 120 K 


3230 ± 130 K 


^phot 


24.2 ± 11.2 pc 


22.7 ±6.8 pc 



a Equin ox J2000.0, epoch 2000 0. 
b From ILepine fc Sharal (|2005l) . These 
authors measure only the combined 
proper motion of the pair, and not indi- 
vidual proper motions for the two stars. 
We have been unable to locate individual 
estimates in the literature. 
c We assume a ±100 K systematic uncer- 
tainty in the effective temperatures. 



only. 

As a check that these objects are indeed nearby dwarfs, 
we show in Figure [2] the position of our source on a V 
versus V — J reduced proper motion ( hereafter, RPM) 
diagra m, reproduced from Fig. 30 of ILepine fc Sharal 
( 2005). Only NLTT 41136 is shown since these authors 
do not provide estimates of magnitudes for NLTT 41135. 
It is not clear if the source was resolved in the surveys 
used to gather their photometry (USNO-B1.0 for V and 
2MASS for J), so we indicate the 0.4 mag uncertainty 
resulting from the possible inclusion of NLTT 41135 in 
either, or both, bands by the parallelogram shaped box 
on the diagram. Nevertheless, it appears the position 
of NLTT 41136 is consistent with that expected for a 
nearby dwarf in RPM. 

We derive photometric estimates of the spectral types 
for both stars using th e average SDSS colors for M-dwarfs 
of iWest et"aT1 ([20051 ). interpolating between the mea- 
sured r — z values in their Table 1. The derived spectral 
types and their uncertainties are reported in our Table 
[4j Note that these uncertainties are largely systematic 
in nature, i.e. the error in relative spectral type between 
the two stars is smaller than the errors we report. We 
have also converted the observed s pectral types to effec- 
tive te mperature using the scale of iKenvon fc Hartmannl 
(1995). We assume a systematic uncert ainty of ±100 K 
in the effective temperatures following iTorres fc Ribasl 
(|200l . 

Finally, photometric estimates of the distance to each 
star were derived using the absolute magnitude relations 
from IWest et"aT1 (|2005[ ) , and comparing to the measured 
magnitudes. These distances differ by < ler, which, in 
combination with their apparent common proper motion, 
strongly supports our conclusion that the two stars are 
physically associated. 




Fig. 2. — V versus V — J r educed proper m o tion diagram for 
LSPM stars, reproduced from Lcpinc & Shara (2005), but plot- 
ted as greyscales and showing the position of NLTT 41136 by the 
parallelogram shaped box. See the text for a discussion of the 
uncertainties in determining the correct position for this object. 

2.8. Ephemeris determination 

The small number of radial velocity measurements 
gathered does not in practice place useful constraints on 
the ephemeris. We therefore derived this entirely from 
photometry, using the discovery data (after applying the 
filtering detailed in £|2.1[) and the times reported in §2] for 
the two well-sampled eclipses. These were modeled with 
a linear ephemeris, fitting the KeplerCam light curve 
model from |j4] simultaneously to the discovery data and 
the two eclipse timings. 

We estimate the uncertainties in the ephemeris pa- 
rameters using the "residual perm utation" bootstrapping 
method (e.g. I Winn et"aT1 12009) to account for the ef- 
fects of the correlated noise (systematics) present in the 
MEarth data on the eclipse times. The light curve was 
first fit adjusting three parameters: the magnitude zero- 
point, ephemeris zero-point io and period P, to compute 
best-fitting values for these parameters. We then pro- 
ceeded to fit the model to new light curves generated by 
taking the original best-fit and adding in cyclic permuta- 
tions of the residuals from that fit. The final parameter 
values and their associated uncertainties were estimated 
using the median and central 68.3% confidence intervals 
of the distributions of the new parameters from these fits, 
over all N — 1 cyclic permutations of the data (where 
N = 993 is the number of data points in the MEarth 
light curve) . In this way, the correlations in the noise are 
taken into account in the results. 

We find that the times of mid-primary eclipse are 
given by tk = to + kP where k is an integer, with 
t = 2455002. 80232±0.00022 (HJD) and P = 2.889475± 
0.000025 days. 

3. RADIAL VELOCITY ANALYSIS 

We fit a single Keplerian orbit model to the observed 
radial velocities using a Levenberg-Marquardt non-linear 
least-squares method. Since the M-dwarf host star 
NLTT 41135 emits most of its light in the infrared and 
the spectral coverage of FIES ends at approximately 7400 
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Fig. 3. — Upper panel: Radial velocity measurements from the 
Nordic Optical Telescope as a function of orbital phase (measuring 
from at the primary eclipse and in units of the orbital period; 
for the assumed circular orbit the secondary eclipse is at 0.5) with 
the best orbital fit overplotted. The 7 velocity of the system has 
been subtracted and the fit assumes a circular orbit, fixing the 
ephemeris to that found by the photometry. The internal error 
bars are plotted, but are not visible due to their small size. Lower 
panel: Phased residuals of the velocities after subtracting the best 
fit model. The rms variation of the residuals is 39 m s" 1 . Note 
that the scale of the lower panel has been expanded to better show 
the residuals. 

A, we were able to use only the reddest orders that were 
not affected by telluric lines. Typically 10 orders were 
used in the analysis. The average internal error for the 
30 minute exposures was 37 m s . We have quadrati- 
cally added a "jitter" term of 17 m s _1 to the velocities 
to yield a reduced x 2 = 1.0. We therefore estimate that 
the true error of the velocities is ~ 41 m s _1 . 

We fixed the mid-eclipse time and period to the values 
found from the photometric analysis of the light curve 
and assumed a circular orbit, which leaves two free pa- 
rameters, namely the velocity amplitude and the 7 ve- 
locity of the system. The radial velocities and resulting 
best-fitting model are shown in Figure [3j and the param- 
eters derived therefrom are presented in Table [5] The 
rms variation of the residuals from the best fit model 
was 39 m s _1 . 

In order to place an upper limit on the eccentricity, 
we also carried out fits where the eccentricity and ar- 
gument of periastron were allowed to vary (whilst en- 
forcing consistency with the photometric ephemeris). 
This yielded an eccentricity of e = 0.007 ± 0.005 and 
K = 13.269 ± 0.062 km s _1 . We can thus constrain the 
eccentricity to be e < 0.02 with 99% confidence. 

4. LIGHT CURVE ANALYSIS 

For detached eclipsing binaries on circular orbits, the 
radial velocity (RV) and light curve models are largely in- 
dependent. Since the radial velocity measurements show 
no evidence for non-zero eccentricity, we assume a circu- 
lar orbit for the purposes of analyzing the light curves. 

Since we do not detect secondary eclipses, and the light 



TABLE 5 

Parameters for the radial velocity model of 
NLTT 41135. 



Parameter 


Value 


7 (km s _1 ) 


41.1 ± 1.2 


Ki (km s" 1 ) 


13.189 ±0.019 


e 


< 0.02 (99% confidence) 



from the secondary is expected to be negligible in the 
optical, with a predicted z-band luminosity ratio of at 
most ~ 10~ 3 between NLTT 41135B and NLTT 41 135 A 
for ages appropriate to an old, field system (see S}6]), we 
assume a dark secondary for the purposes of modeling 
the light curve. This allows the use of simple models 
developed for fitting transiting exoplane t light curves. 

We use the analytic transit curves of lMandelfc' Agol 
(2002) to model the KeplerCam primary eclipse light 
curve, modeling limb darkening with a standard 
quad ratic law (e.g. iClaret fc Gimenezll 19901 Ivan Hammd 
1993) of the form: 

|M = 1- Ul (l- M )- U2 (l- M )2 (1) 

where Ix(fJ-) expresses the monochromatic intensity at 
wavelength A for a point on the surface at angle 9 from 
the normal, /1 = cos 9, and u\ and ui are free parameters. 

We accounted for the dilution of the eclipse depth by 
the presence of NLTT 41136 in the photometric apertures 
using the measured z-band light ratio (L3) from ^2.21 
Since the light ratio and the light curves themselves were 
both observed through z filters, the error introduced by 
doing this should be small, and we account for its effect in 
the final uncertainties by assuming a conservative error 
on L3. 

After accounting for the known system ephemeris, the 
model has five free parameters: the ratio of component 
radii R2/R1, orbital inclination i, the semimajor axis 
divided by the primary radius, a/Ri, and two limb dark- 
ening parameters Ui, u%. In practice, we introduce three 
additional parameters: the normalization (out of eclipse 
magnitude) zq, a linear term in airmass k(X — 1) (where 
X is airmass) to account for any residual differential at- 
mospheric extinction effects, and we also allow for a tim- 
ing offset At to account for any error in the ephemeris 
or deviations in the eclipse time due to orbital perturba- 
tions. 

While it is usually possible (with sufficiently high- 
quality data) to fit for all three geometric parameters 
(R2/R1, i and a/R\) of a total eclipse, it is important to 
note that this is not possible for a grazing eclipse, because 
the second and third contact points are no longer seen 
in the light curve. In the simple case with no limb dark- 
ening and a dark secondary, this can be understood as 
follows. For a total eclipse, it is straightforward to show 
that the eclipse depth is determined by the area ratio 
between the two objects, and is thus equal to (R2/R1) 2 ■ 
The eclipse duration, and the duration (or equivalently, 
slope) of the partial phases (first to second contact, and 
third to f ourth contact) are determined by i and a/Ri 
(e.g. see Sea ger fc Mallen-Ornelasl [2003 for a detailed 
derivation). In the grazing case, although we still know 
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the eclipse duration and the slope of the partial phases, 
we do not know what the level of the "bottom" of a total 
eclipse would be. Therefore, the parameter R2/R1 is es- 
sentially undetermined without external information on 
the other two. 

Since we are primarily interested in i and a/R\, we 
adopt a weak Bayesian prior on the most poorly deter- 
mined parameter, R2/R1, and marginalize over this pa- 
rameter to determine the distribution of possible values 
of the other two. In addition, the grazing light curve only 
constrains the limb darkening parameters very weakly, so 
we fix these a t values appropriate for the z passband from 
IClaretl ([2004) using the effective temperature from Table 
[4]and assuming log g = 5.0 and solar metallicity. The 
theoretical uncertainty in limb darkening parameters for 
the z passband should be negligible compared to the in- 
trinsic degeneracies in the modeling of the grazing eclipse 
for the present case. 

In order to derive parameters and reliable error es- 
timates, we adopt a variant of the popular Markov 
Chain Monte Carlo analysis f requently applied for tran- 
siting exoplanet s ystems (e.g. [Tegmark 2 0041 : iFordl [20051 : 
iWinn et al.ll2007l ) . The majority of thes e analyses use the 
standard Metropoli s-Hastings method ([Metropolis et all 
119531 : [Hastings ll970D with Gibbs sampler, which has the 
advantage of being extremely simple to implement. We 
briefly describe this method here (as typically used for 
parameter estimation in astronomy), and then summa- 
rize the modifications made in our implementation. 

Starting from an initial point in parameter space, the 
Metropolis-Hastings algorithm takes the most recent set 
of parameters and perturbs one or more parameters by 
a random Gaussian deviate. Gibbs sampling perturbs 
one parameter at a time, cycling around the parameters. 
If the perturbed parameter set has a lower \ 2 than its 
progenitor, it is accepted as a new point in the chain. 
If it has a larger \ 2 i it is accepted with a probability 
exp(— Ax 2 /2). If it is not accepted, the original point is 
repeated in the chain. The size of the perturbations are 
usually adjusted by manual iteration so that 20 — 30% of 
the proposed points are accepted. 

We adopt the Adaptive Metropolis algorithm of 
lHaario et al.l ((2001) , which dynamically updates a Gaus- 
sian proposal distribution during the simulation using 
the empirical covariance matrix, thus largely eliminating 
the need for initially tuning the proposal distributions to 
obtain the correct acceptance rates. This is by defini- 
tion non-Markovian, but the updates are done in a way 
that has been shown to maintain the correct ergodicity 
properties of the chain. 

We used a standard Levenberg-Marquardt algorithm^] 
to provide the initial parameter and covariance estimates 
to start the chain, and used the resulting fit to re-scale 
the observational errors such that the reduced \ 2 was 
equal to unity. 

For parameter estimation, we used a chain of 10 7 
points, discarding the first 20% of these in order to ensure 
the chain had converged. The correlation lengths for all 
parameters were < 100 points. We enforced the observed 
luminosity ratio between NLTT 41136 and NLTT 41135, 
of L 3 = 2.1 ± 0.2 as discussed in 

As discussed earlier in this section, we have assumed 

4 http: //www. ics . forth. gr/~lourakis/levmar/ 



a prior on the radius ratio, to break the degeneracy be- 
tween the parameters in the light curve model. With- 
out any prior, the allowed values of R2/R1 are bounded 
only below (corresponding to the limiting case where the 
eclipse switches from being grazing to total), and are 
essentially unbounded above despite R2/R1 > 1 being 
physically extremely unlikely. 

From the radial velocities, we are left with little doubt 
that t he secondary lies below th e hydrogen burning limit 
(e.g. [Chab rier fc Baraffe) 120001 and references therein) 
and is therefore a brown dwarf. The physics of this tran- 
sition is quite well understood, and except for deuterium 
burning, the only source of energy for these objects is 
gravitational contraction. Theoretical models predict 
that all objects with masses i n this range should have 
radii < 0.11 R© after 300 Myr ([Burrows et al.lll997[ ). im- 
plying r R 2 /E J _^0.6^The measured radii of OGLE-TR- 
122b ([Pont et al.l 12001 and CoRoT-3b ([Deleuil et al.1 
2008) also corroborate this theoretical argument. We can 
place an analogous lower limit, for example, secondary 
sizes smaller than the planet Neptune would require im- 
plausibly high densities at the measured secondary mass, 
implying R 2 /R\ > 0.17. 

We therefore adopt a Gaussian prior, R2/R1 — 0.4 ± 
0.3, the value of which was derived by taking the pre- 
dicted radius ratio from the stellar and brown dwarf 
models given the estimated masses of the two compo- 
nents, but also encapsulates the argument made in the 
previous paragraph. We assume a conservative error to 
account for the uncertainty in this estimate. An isotropic 
prior (i.e. uniform in cost) was assumed on the orbital 
inclination, and uniform priors were assumed on zq, k 
and R\/a. 

Our final parameter values and their associated un- 
certainties were estimated using the median and central 
68.3% confidence intervals of the marginalized posterior 
probability distributions for each parameter, and are re- 
ported in Table [U Figure 0] shows the best- fitting model 
overplotted on the light curves, and in Figure [5] we show 
representative two-dimensional distributions of the pa- 
rameters of the model to illustrate the discussion of the 
degeneracies in the light curve modeling from this sec- 
tion. It is clear from the table that despite the inher- 
ent degeneracies in the light curve model, we do indeed 
obtain acceptable bounds on the inclination and a/R\, 
especially recalling that the inclination appears as sinz 
when we come to interpret the radial velocity measure- 
ments. 

Note that the choice of a linear term in airmass k(X — 1) 
to represent the variations in the out-of-eclipse baseline is 
somewhat arbitrary. While it is reasonable to expect at- 
mospheric extinction to affect the differential photometry 
at the observed level (especially given the much redder 
color of our target relative to the comparison stars), it 
is equally possible the observed "slope" is due to stellar 
variability. This means in practice that the depth in- 
formation used in the model fits is uncertain due to the 
unknown contribution of stellar spots. 

Assuming a sinusoidal out of eclipse modulation with 
the same period as the orbit, the measured value of k 
would correspond to a semi-amplitude of > 0.01 mag, or 
correspondingly, > 0.005 mag if the period was half an 
orbit (e.g. ellipsoidal variation) . While the MEarth pho- 
tometry places limits on the actual level of out of eclipse 
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TABLE 6 

Parameters for the light curve model of 
NLTT 41135. 



Parameter Value 



to (HJD) 2455002.80232 ± 0.00022 

P (days) 2.889475 ± 0.000025 

u± 0.0169 (assumed) 

U2 0.6976 (assumed) 

L 3 2.10 ±0.20 (assumed) 

At (HJD) 0.00006 ± 0.00023 

k (mag/airmass) -0.0062 ± 0.0023 

i (deg) 87.42±g;" 

a/Ri 24.60tjgg 

Ba/fli 0.48«;^ 
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Fig. 4. — Top panel: phase-folded 2-band KeplerCam light curve 
with our best-fitting light curve model over-plotted. Center panel: 
as above, with the linear airmass term subtracted out before plot- 
ting. Bottom panel: residual after subtracting the model from the 
data. 

variation (see i)2.ip . and thus on the possible phases 
and periods corresponding to the measured fc, the un- 
explained out of eclipse variations seen there may indeed 
be consistent with the "slope" seen in the KeplerCam 
data. 

We use the resolved light curve from §2.2! to verify the 
parameters we have derived from the KeplerCam curve, 
shown in Figure [6l We re-fit for the three parameters zq, 
k and At, and assumed L% — (since the light curve was 
resolved), but all other parameters were held fixed from 
the KeplerCam analysis. The fit was carried out using 
3cr clipping, in order to reduce the effect of the large 
number of photometric outliers on the fit results, and we 
re-scaled the observational errors such that the reduced 
X 2 of the out of eclipse parts of the light curve was equal 
to unity. The reduced x 2 01 the mu fit (including the 
eclipse) was 1.04. 

As an additional check, we have also performed a 
simultaneous fit of the KeplerCam and OPTIC light 
curves. Separate normalization, k and At values were 
allowed for each of the two curves, and we fit for a single 
set of geometric parameters (k, i, and a/R\) and for L3 
(removing the prior used in the previous analyses and 
instead fitting for this parameter). We used the same 



TABLE 7 

Measured mid-eclipse times for NLTT 41135. 



HJD a 


(O - C) (s) Cycle b 


Instrument 


2454999.91269 


-13 ±37 -1 


OPTIC 


2455002.80238 


5 ±20 


KeplerCam 



a Heliocentric Julian Date of mid-eclipse. All 
HJD values reported in this paper are in the 
UTC time-system. 

b Cycle number, counting from at the pri- 
mary eclipse at to, in units of the orbital pe- 
riod. 

clipping and rescaling of errors carried over from the sep- 
arate analyses, which should correct for any differences 
in the under-estimation of the observational uncertainties 
between the two curves and produce the correct relative 
weights. 

The values of i, R2/R1 and a/Ri from this analysis 
differ by <C la from those given in Table O providing 
further verification of the results. We derive L3 = 2.15 ± 
0.16, which is consistent with the value assumed in the 
analysis of the KeplerCam light curve to < la. The 99% 
confidence interval from this fit is 1.80 < L3 < 2.52. 

Finally, we report the mid-eclipse times from our anal- 
ysis of the two well-sampled eclipses from the unresolved 
KeplerCam and the resolved OPTIC light curves in Table 

m 

We now proceed in the next section to interpret the 
parameters derived from the radial velocity and light 
curve modeling in terms of the physical parameters 
of the M-dwarf NLTT 41 135 A and the brown dwarf 
NLTT 41135B. 

5. PROPERTIES OF NLTT 41135A AND SYSTEM 
PHYSICAL PARAMETERS 

Using the 7 velocity from Table [5] and the mea- 
sured positions, proper motions and photometric par- 
allax from Table SJ we estimate (U, V, W) = (-53.8 ± 
6.2, -5.0 ± 3.0, 4.9 ± 5 . 9) km s" 1 using the method of 
Uohnson fc Soderbloml (|1987l) . where we adopt the defi- 
nition that positive U values are away from the Galactic 
center. This places NLT T 41135 in the old Galactic disk 
population as defined by Leggctt (1992). While kinemat- 
ics provide relatively weak constraints on the stellar ages, 
these are typically a few to 10 Gyr for old disk stars, and 
it is thus highly likely that NLTT 41135 is older than a 
few Gyr. 

In order to derive an estimate of the m ass of 
NLTT 41135A, we adopt the methodology of iTorres! 
(|2007f ). although the available information in our case 
is more limited since we do not have a reliable measure- 
ment of the parallax or near-infrared magnitudes. 

The quantity a/R\ from the light curve model in iJH 
provides a measurement of the mean density of the pri- 
mary, p\ using Newton's version of Kepler's third law: 

/ a \ 3 _ G(Mi + M 2 ) G(l + g) 47rpj 
\rJ i??r2 2 VP 3 U 

where f2 = 2ir/P is the angular frequency of the or- 
bit. We explicitly note that the secondary mass is not 
negligible in the present case. 
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Fig. 5. — Two-dimensional grcyscale histograms of the posterior probability density functions for selected parameter pairs from our Monte 
Carlo analysis. The mapping from probability density to the intensity of the greyscales was a square root function, to compress the dynamic 
range and allow the tails of the distribution to be more readily distinguished. 





Fig. 6. — Phase-folded resolved z-band light curve from the OP- 
TIC data with our best-fitting light curve model over-plotted. Pan- 
els as Figure [4] Points that were excluded from the fit by our clip- 
ping procedure are colored red in the online version of the figure. 



We use this estimate of the stellar density in conjunc- 
tion with stellar evolution m odels to derive M\. F or this 
work, we adopt the tracks of Baraffc et al.l (|1998f) . 

It is well- established from stud i es of eclipsing binary 
stars (e.g. lLopez Morales! 120071: iChabrier etall 120071 1 
that stellar evolutionary models tend to underpredict the 
radii of M-dwarfs at a given mass. The bolometric lumi- 
nosities from the same models do not seem to be seriously 
in error (e.g . iDelfosse et al.|[2000t rTbrres fc Ribasl [20021: 
iTorres et all 120061: ITorresI I2007D . We therefore follow 



iTorresI (|2007j ) and introduce a multiplicative correction 
factor /3 to the radii, while at the same time applying a 
factor to the effective temperatures to conserve the 

bolometric lumin osity. Studies of eclipsing binaries (e.g. 
iRibas et al.ll2008D and the transiting planet host GJ 436 
( Torresll2007t ) indicate that typical values of (3 ~ 1.1. 

In practice, we fit for the observed stellar density and 
effective temperature (from Tabled]), given a fixed metal- 
licity, and allowing the age to vary. We do not use the 
observed optical colors as an additional constraint, be- 
cause these are known to be po orly- reproduced by the 
models (e.g. iBaraffe et al.l ll99g). This is thought to be 
the result of a missing source of molecular opacity in the 
NextGen model atmospheres. 

The kinematic information presented in this section 
indicates that our target is an old, field M-dwarf. Wc 
therefore explore a range of ages from 1 — 10 Gyr, not- 
ing that this has little effect on the derived parameters 
in practice, because M-dwarfs in this age range evolve 
relatively slowly. The metallicity is also unknown, but 
typical values for old disk stars are [Fe/H] ~ —0.5 (e.g. 
Leggett 1992). We derive parameters for both this value 
and for solar metallicity, to give an idea of the likely range 
resulting from this source of uncertainty, which does have 
a significant effect on the stellar parameters. 

Our derived physical parameters for NLTT 41135A, 
the orbital parameters, and the inferred mass of 
NLTT 41135B, are given in Table M The value of 
/3 indicates essentially no inflation of the radius of 
NLTT 41 135 A relative to the models, within the ob- 
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TABLE 8 

Derived physical parameters for the 
NLTT 41135 system. 



Parameter Value 

[Fe/H] = 0.0 [Fe/H] = -0.5 



Mi (M ) 




164+ - 021 


Ri (R©) 


210+ ' 016 
u.ziu_ 014 


201+ - 014 
u.zui_ 013 


Li (L ) 


0043+ 0012 
u.uuia_ 0010 


0039+ 0012 
U - UUOJ -0.0009 


log 9i 


5 062+ 033 

O.UDZ_ 034 


c n79 +0.033 
°- u,z -0.034 





1 021+ ' 057 
i.uzi_ 053 


1 ova+o.osf) 

i-U'»_0.055 


a (Rq) 


°' io -0.20 


4Q3+019 


1 


n 171 +0.008 
u - 1 1 -0.008 


n ion+0-008 


M 2 (M Jup ) 


oq y+2.8 
°°- 1 -2.6 


30.9± 2 ; 4 



servational errors, but this is not surprising since the 
propagated error in the effective temperature measure- 
ment yields only a weak constraint on this parameter, 
with the value of 1.1 found in other studies also being 
reasonably consistent with the measurements. 

Given the limited constraints available on the proper- 
ties of NLTT 41135A, it is perhaps not surprising that 
the knowledge of the mass of the brown dwarf secondary, 
NLTT 41135B, is presently limited by the uncertainty in 
the primary mass. One of the most reliable methods to 
constrain single M-dwarf masses in this range is to use 
the mass-lumin osity relation, either empirical determi- 
nations such as Dclfoss e~et al.l ([2000) or the models. We 
therefore suggest that the most profitable way to pro- 
ceed would be to obtain an accurate trigonometric paral- 
lax, and resolved, near-infrared apparent magnitudes for 
NLTT 41135 and NLTT 41136. Of the available pass- 
bands, K woul d appe ar to be the logical choice, since the 
iDelfosse et al.l (|2000D K-band relation is one of the best- 
constrained and shows the smallest scatter. This should 
allow the uncertainty in the primary mass to be reduced 
to ~ 10%, and permit a test of the models of M-dwarfs 
in addition to improving the secondary parameters. We 
also note that the passband-integrated luminosities of 
M-dwarfs can be used to estimate the metallicity (e.g. 
I Johnson k Appll2009l ). 

Since NLTT 41135 is a member of a presumed visual 
binary system, ultimately it may be possible to mea- 
sure its dynamical mass from the orbit of the visual bi- 
nary. This would allow a model-independent radius to 
be derived for NLTT 41 135 A, permitting a direct test 
of M-dwarf evolution models, and would also provide a 
model-independent dynamical mass for NLTT 41135B. 
We discuss this possibility in §T\ 

6. PROSPECTS FOR CONSTRAINING BROWN DWARF 
EVOLUTION MODELS 

While we estimate the mass of NLTT 41135B to 
10 — 15% (limited by the uncertain metallicity and accu- 
racy of determination of the mass of the M-dwarf host), 
the lack of a secondary eclipse and the detection of only 
one object in the spectrum mean that this is the only 
parameter we can presently determine. 

The most straightforward approach to obtaining a sec- 
ond fundamental parameter for the brown dwarf is to 
attempt to detect a secondary eclipse. The ratio of the 



TABLE 9 

Predicted luminosity ratios, surface brightness ratios, 
and secondary eclipse depths for nltt 41135 with 

Ag = 0.0. 



Age T 2 Band L 2 /Li J d 2 a 

(K) (%) 



1 Gyr 1067 z 8.60 x 10" 4 3.94 x 10" 3 0.03 

J 3.51 x 10~ 3 1.61 x 10" 2 0.12 

H 1.97 x 10~ 3 9.04 x 10~ 3 0.068 

K 1.59 x 10~ 3 7.29 x 10~ 3 0.055 

L' 1.04 x 10~ 2 4.76 x 10~ 2 0.36 

M 1.73 x 10~ 2 7.95 x 10~ 2 0.6 

5 Gyr 643 z 1.17 x KT 4 6.24 x 10~ 4 0.0047 

J 1.84 x 10" 4 9.81 x 10~ 4 0.0074 

H 1.09 x 10~ 4 5.82 x 10~ 4 0.0044 

K 2.04 x 10" 5 1.09 x 10~ 4 0.00082 

L' 1.18 x 10" 3 6.31 x 10" 3 0.047 

M 4.26 x 10~ 3 2.27 x 10~ 2 0.17 

10 Gyr 533 z 4.64 x lO -5 2.61 x 10" 4 0.002 

J 5.19 x 10~ 5 2.92 x 10" 4 0.0022 

H 3.14 x 10~ 5 1.77 x 10" 4 0.0013 

K 2.15 x 10~ 6 1.21 x 10~ 5 9.1 x 10" 5 

L' 4.76 x 10~ 4 2.68 x 10~ 3 0.02 

M 2.31 x 10~ 3 1.30 x 10~ 2 0.097 



a Predicted secondary eclipse depth. We have assumed 
a primary eclipse depth of 7.5% based on the light curves 
shown in SjD 

primary to secondary eclipse depths in a given pass- 
band is approximately the ratio of central surface bright- 
nesses of the two objects, J = L 2 i?i/Lii?|, where all 
these quantities are measured for the specific passband 
in question. The division by the squared radius ratio 
makes this a more favorable prospect than direct detec- 
tion of light from the secondary e.g. in a spectrum, and 
many successful detections of shallow secondary eclipses 
have been achieved for exoplanet systems, especially in 
the mid-infrared using the Spitzer space telescope (e.g. 
iCharbonneau et al.ll2005l: iDeming et all [20051 [20061) . 

Table [HI gives predicted luminosity ratios, surface 
brightness ratios, and secondary eclipse depths assum- 
ing the value of Mi from Table [5] and using brown dwarf 
evolution models to predict the surface brightness of 
NLTT 41135B i n near-infrared band passes. We use the 
COND models of iBaraffe et al.1 ([2003) for the secondary, 
and the NextGen models discussed in |J5] for the pri- 
mary. We do not compute predictions for bluer band- 
passes than z due to known problems w ith both sets 
of m odels, as mentioned earlier (see also IBaraffe et al.l 
120031 for discussion of problems with the I-band mag- 
nitudes in the COND models), but the predicted eclipse 
depths are even shallower at these wavelengths due to 
the rapidly declining flux from the brown dwarf, so de- 
tection here seems unlikely. Three representative ages, 1, 
5 and 10 Gyr, are shown in the table, with the kinematic 
evidence discussed in fj5] favoring the old end of this age 
range. 

For the purposes of these calculations, we have ne- 
glected stellar insolation. The zero albedo equilibrium 
temperature for NLTT 41135B is approximately 460 K, 
so it is likely that we underestimate the surface bright- 
nesses, especially at the two oldest ages where insola- 
tion will represent a significant contribution to the brown 
dwarf energy budget. 
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TABLE 10 

Predicted luminosity ratios, surface brightness 
ratios, and secondary eclipse depths for nltt 41135 

WITH A g = 1.0. 





T> 


Band 


L 2 /L 1 




./ 




d,2 a 




(K) 










(%) 


1 Gyr 


1067 


z 


1.22 X 10" 


-3 


5.59 x 10- 


-3 


0.042 






J 


3.87 x 10" 


-3 


1.78 x 10" 


-2 


0.13 






H 


2.33 x 10" 


-3 


1.07 x 10" 


-2 


0.08 






K 


1.95 x 10" 


-3 


8.95 x 10- 


-3 


0.067 






V 


1.07 x 10" 


-2 


4.92 x 10- 


-2 


0.37 






M 


1.77 x 10" 


-2 


8.11 x 10- 


-2 


0.61 


5 Gyr 


643 


z 


4.27 x 10" 


-4 


2.28 x 10- 


■3 


0.017 






J 


4.94 x 10" 


-4 


2.63 x 10- 


3 


0.02 






H 


4.19 x 10" 


-4 


2.23 x 10- 


3 


0.017 






K 


3.30 x 10" 


-4 


1.76 x 10- 


3 


0.013 






L' 


1.49 x 10" 


-3 


7.96 x 10- 


3 


0.06 






M 


4.57 x 10" 


-3 


2.44 x 10- 


2 


0.18 


10 Gyr 


533 


z 


3.40 x 10" 


-4 


1.91 x 10" 


•A 


0.014 






J 


3.45 x 10" 


-4 


1.94 x 10- 


3 


0.015 






H 


3.25 x 10" 


-4 


1.83 x 10- 


-3 


0.014 






K 


2.96 x 10" 


-4 


1.66 x 10- 


3 


0.012 






V 


7.69 x 10" 


-4 


4.33 x 10- 


-3 


0.032 






M 


2.60 x 10" 


-3 


1.46 x 10- 


2 


0.11 



a Predicted secondary eclipse depth. We have as- 
sumed a primary eclipse depth of 7.5% based on the 
light curves shown in fJU 

It is also important to consider the validity of the as- 
sumption of no "reflected light" from the brown dwarf, 
since although its temperature is predicted to be very 
low, the radius should be only ~ 1/2 — 1/3 that of the 
primary. For brown dwarfs with effective temperatures 
< 1300 K, complete grav itational sett l ing o f any dust 
grains should occur (e.g. lAllard et~ al. 2001), which is 
the physics assumed by the COND models we have used, 
and should give low albedo, since the atmosphere consists 
mostly of absorbing species. We compute an additional 
set of predictions (shown in Table I10[) assuming an (un- 
realistically large) value for the geometric albedo (A g ) 
equal to unity to illustrate an maximum size of this ef- 
fect, which are reproduced in Table [TO] Comparing the 
results in the two tables indicates that reflection should 
have only a very minor effect on our predictions for the 
bandpasses where the secondary eclipse is most likely to 
be observable. 

It is clear from the tables that the best prospects for 
detection of the secondary eclipse lie at the longest wave- 
lengths (L/ and M), and will likely require space-based 
observations in order to obtain adequate precision. J 
and to some extent H show more favorable eclipse depths 
than K due to the familiar blueward shift of the J — K 
and H — K colors for T-dwarfs resulting from H2 and 
CH4 absorption features. 

Assuming M-band to be a reasonable proxy for the 
longest-wavelength (4.5 fxm) Warm Spitzer channel, this 
eclipse might be detectable, although it is likely that 
Spitzer's poor angular resolution would result in the im- 
ages of NLTT 41135 and 41136 being blended, thus di- 
luting the eclipse depths by a factor of ~ 3. Whilst the 

5 Upper limits from observations of secondary eclipses and phase 
curves in the optical for Hot Jupiter exoplanets, which are thought 
to have similar atmospheric const ituents, indicate that A g < 0.1 in 
the visible (e.g. IRowe et ani2008l 1. 



eclipses should still be detectable for the youngest ages 
we have considered, this might become extremely chal- 
lenging at the more likely age of ~ 10 Gyr. Due to the 
rapid decline in J-band as a function of age, a similar 
statement applies for the prospects of detection using 
near-infrared instruments aboard the Hubble Space Tele- 
scope. 

With the James Webb Space Telescope, the secondary 
eclipse detection should be straightforward, and the pre- 
dicted angular resolution is sufficient to resolve the M- 
dwarf pair even out to the mid-infrared. 

It is interesting to speculate that if the secondary 
eclipse was total (which we stress is very unlikely given 
the essentially circular orbit and grazing geometry of 
the primary eclipse), this would yield the bandpass- 
integrated luminosity of the brown dwarf directly from 
the eclipse depth, and would break the degeneracy inher- 
ent in the modeling of the grazing eclipses in conjunction 
with the radial velocity, yielding also the radii of both 
components. 

Finally, we note that it may be possible to derive a 
relatively crude constraint on the radius of the brown 
dwarf even in the most likely case where both eclipses 
are grazing. Recall from fj4] that we assumed a prior on 
the system radius ratio in order to determine a/ R\. How- 
ever if we knew the radius of the M-dwarf, a/R\ would 
be already completely determined and could be fixed 
in the model, allowing us to instead constrain R2/R\. 
This could be done by measuring the angular diameter 
of NLTT 41135, in conjunction with a parallax and a 
K-band magnitude measurement to determine its mass. 
Due to the remaining strong degeneracies between R2 / R\ 
and i it is unlikely that this measurement would yield 
R2/R1 to any better than ±0.1 (see Figure [5]). 

7. THE NLTT 41135 / 41136 VISUAL BINARY SYSTEM 

Using the iBaraffe et all (|1998t ) models, we estimate a 
mass ratio of ~ 0.8 for the M-dwarf visual pair from the 
effective temperatures in Table 0] (assuming [Fe/H] = 
—0.5, 5 Gyr age and the value of (3 from Table [8]). Assum- 
ing the mass of NLTT 41135 from Table [U this implies 
the mass of NLTT 41136 is approximately 0.21 M Q , and 
the total system mass is therefore M tot ~ 0.37 M Q . The 
angular separation of 2'.'40 corresponds to a projected 
physical separation of 55 AU at 23 pc distance. Given 
this information, we can estimate a minimum orbital pe- 
riod of 680 years using Newton's version of Kepler's third 
law. 

Given the long orbital period, it will be challenging to 
measure the parameters of the astrometric orbit unless 
its orientation or eccentricity are extremely favorable, al- 
though the motion of the components relative to their 
presumed common proper motion could be as large as a 
few tens of milliarcseconds per year especially if the or- 
bit is viewed face-on. The maximum radial velocity dif- 
ference between the pair corresponding to this period is 
2.4 km s _1 (assuming an edge-on orbit), which is reason- 
ably consistent with our measurements from the TRES 
spectroscopy (see £|2.5[) provided the measurements were 
taken at less than the maximum separation or the orbit 
is not edge-on. 

8. DISCUSSION 
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We have reported the discovery of a hierarchical triple 
system containing two M-dwarfs and a brown dwarf or- 
biting the less massive of the pair in an eclipsing system. 
The masses of the components are approximately in the 
ratio 8:6:1. 

The existence of this object poses a challenge to theo- 
ries of brown dwarf formation. The configuration of the 
system strongly resembles a scaled-down (in mass) ver- 
sion of a relatively common configuration for solar-type 
stars. The mass ratio of only 6 : 1 for the eclipsing binary 
component of the system would seem to be difficult to 
produce by disk fragmentation, because it would require 
an extraordinarily massive disk surrounding the forming 
M-dwarf core. 

Given this difficulty it seems likely that this object 
may have formed in the same manner as its higher-mass 
analogs: directly from gravitational collapse of an over- 
density in a molecular cloud. It is interesting to postulate 
that if many such systems are made by the star forma- 
tion process, there will be some that are dynamically un- 
stable, and may eject one or more of their components. 
As propose d in the ejection hypothe sis of brown dwarf 
formation (jReipurth fc Clarka 12001). this may offer a 
mechanism to produce single brown dwarfs, although as 
discussed in SJTJ evidence such as the existence of single 
brown dwarfs with discs argues against this being the 
only mechanism for brown dwarf formation. 

It is interesting to compare the secondary of the present 
system with the secondary of the CoRo T-3 system, a 
22 M .Tnp brown dwarf orbiting an F3V star (jDeleuil et al.1 
2008) . The brown dwarfs in the two systems have similar 
masses, and their orbital periods differ by only 50%, yet 
CoRoT-3b orbits a star of luminosity 10 3 times greater 
than does NLTT 41135B, and thus experiences a 500 
times greater (bolometric) stellar insolation. This means 
that while the energy budget of NLTT 41135B is most 
likely dominated by internal energy sources, CoRoT- 
3b will be dominated by insolation. It would thus be 
extremely interesting to compare the observed physi- 
cal properties of these objects to obtain insights into 
brown dwarf physics, for example the efficiency of re- 
distribution of energy from the day-side to the night- 
side of C oRoT-3b, as has bee n done for extrasolar plan- 
ets (e.g. iKnutson et al.ll2007l ). We also note that these 
objects may have formed through different mechanisms, 
since the mass ratio for the CoRoT-3 system of q = 0.015 
(an order of magnitude lower than for NLTT 41135B) 
may be difficult to produce via conventional mechanisms 
thought to be responsible for binary star formation. 

NLTT 41135B is also of comparable mass to 
the secondary in the dou ble-lined brown dwarf EB 
2MASS J05352184-0546085 (jStassun et all 120061 [2001 . 
although it is of much older age and orbits a main se- 
quence star. Nonetheless, these objects are probably 
quite close to being evolutionary analogues of one an- 
other, and both have an energy budget most likely dom- 
inated by internal energy sources. 

The present system provides a dynamical mass for a 
field brown dwarf, and offers the potential to allow mod- 
els of brown dwarf evolution to be tested at extremely old 
ages, where dynamical constraints are scarce. In order 
to test the theory, we must attempt to determine or con- 
strain the system age, in addition to measuring another 
fundamental property of the brown dwarf. The most ac- 



cessible measurement is the secondary eclipse, which in 
conjunction with the primary eclipse yields a measure- 
ment of the central surface brightness ratio of the two 
components from the ratio of the eclipse depths. 

Finally, we note that a major limitation in the determi- 
nation of the brown dwarf mass and M-dwarf parameters 
is the uncertainty in the mass of NLTT 41135A. We sug- 
gest that the best prospect for improvement would be 
to obtain an accurate trigonometric parallax for this ob- 
ject, in conjunction with resolved near-infrared apparent 
magnitude measurements (preferably in K-band). This 
should yield an estimate for the metallicity, and thus the 
primary mass to ~ 10%, giving an improved solution for 
the eclipsing binary to more tightly constrain models of 
low-mass stars and brown dwarfs. 
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